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We present data from the ALPHA Collaboration about lattice calculation of SU(2) pure-gauge running coupling 
constant, obtained with two different definitions of the coupling itself, which show universality of the continuum 
limit and clarify the applicability of renormalized perturbation theory. 



1. INTRODUCTION 

Computation oiagtrong by means of the finite- 
size technology of ref. [1] have been previously 
performed in the pure SU(2) and SU(3) gauge 
theory [2-7]. Using the massive computer power 
provided by the APE machines at Rome II and 
DESY, we were able to extend these calculations 
in several directions. Here and in a much more 
detailed manner in ref. [8], we discuss such deve- 
lopments. 

2. THE TWO COUPLINGS 

We will not enter into too many details, refer- 
ring the interested reader to [8]. We will only 
recall that the Schrodinger functional involves an 
integration over all gauge fields with spatial pe- 
riodic boundary conditions and fixed boundary 
conditions in time: it may be seen as the propa- 
gation kernel for going from the initial gauge field 
configuration C at time to the final configura- 
tion C" at time L. Taking the derivative of the 
effective action respect to a parameter 77 which en- 
ter in the definitions of classical boundary fields 
C and C", we can define a renormalized coupling 
which we call asF- The only physical scale in the 
system is the box size L, so we can consider thw 
R.G. flow of the coupling arising from a change 
in this scale. 

The other coupling is defined as the ratio of cor- 
relations of Polyakov loops, in a box with twisted- 



periodic boundary conditions. The correlations 
are taken in the time direction at a distance L/2, 
to ensure that the only scale entering in the def- 
inition is once again L. Correlations of Polyakov 
loop winding in one of the two twisted directions 
start in perturbation theory as Qq, while those in 
the periodic directions start as 1. Ratio of such 
correlations, with a normalization prefactor com- 
ing from a tree-level perturbative computation, 
allow us to define axp- 

At this point we would stress the fact that the 
infrared behaviours of asF and axp are com- 
pletely different since the latter diverges expo- 
nentially in the q —f limit, while the former 
is expected to approach a computed constant. 

For both couplings we performed a perturba- 
tive computation giving us the relation between 
the coupling themselves and other more usual def- 
inition: for example 
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asF + 0.9433 x {asF^ 



(1) 



with a second order coefficient which is under 
computation, and 
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utp - 0.5584 X {aTpY 



(2) 



3. SIMULATIONS 
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In our simulations we used a standard hybrid 
over-relaxation algorithm with Nor microcanon- 
ical sweeps followed by one heat-bath sweep. In 
all cases the high-quality random number gener- 
ator of ref. [9] was used. 



Surprisingly enough, we found that the inte- 
grated autocorrelation times of our observables 
depend strongly on the order in which each point 
of the lattice is visited during one sweep. The op- 
timal choice was to update all links in a certain 
direction ji for all lattice points, then switching to 
the next value of ji and so on. The other choice, 
i.e. updating all four positive-directed links of 
a site before passing to the next site, results in a 
slightly more optimized code but in a much higher 
autocorrelation time (a factor 4 for SF and a fac- 
tor 2 for TP, on moderately-sized lattices). 

By studying autocorrelation times and rela- 
tive variances of our observables, we find that 
the overall computational effort scales roughly as 
(L/a)^ when one approaches the continuum limit 
at fixed L for both couplings, with a bigger co- 
efficient (about one order of magnitude) for TP 



4. EXTRAPOLATION TO THE CON- 
TINUUM LIMIT 

The precision that can be reached for the two 
couplings at a given value oi L/a is quite different, 
and large lattices are very expensive for ax p. On 
the other hand lattice artifact are also different in 
the two cases: in particular for axp the leading 
corrections are expected to be 0[[a/ L)^), while 
for asF the presence of "source-walls" at the time 
edges of the lattice allows for terms 0{a/L). 

In TP case we found that not-too-big lattices 
(up to 16'*) are enough for a safe extrapolation to 
the continuum limit, while in SF case we found 
that the use of Symanzik improvementent was of 
some help in reducing 0{a/ L) lattice artifact. 

We consider also a perturbative improvement 
of the observables. For SF case we can build a 
2-loops improved observable 



i:fp{2,u,alL) 



^(2), 



i:sF{2,u,a/L) 



1 + 6i{a/L)u + 62{a/L)u^ 



(3) 



where 11gp{2,u,a/L) is the lattice step scaling 
function with a scale factor s = 2 computed at 
g^p = u, and 6i[a/L) are perturbative correc- 
tions. For TP case we can only build the 1-loop 
improved observable. 



This improvement is indeed very efficient in re- 
ducing the lattice artifacts for S5/, while for Syp 
at the two largest value of the coupling the non- 
perturbative effects start to win over the improve- 
ment, leading to higher values of lattice artifacts. 

5. RENORMALIZED PERTURBATION 
THEORY 

Let us now try a quantitative comparison be- 
tween the two renormalized couplings: we want to 
obtain a non-perturbative universal relation be- 
tween aspiq) and aTp{q). To achieve this goal, 
we decided to fix aspiL) = 0.16535 for several 
values of a/L, and at the corresponding values 
of P we compute aTp{L). Then we can safely 
extrapolate to the continuum limit, keeping only 
the points for L/a > 8 and allowing for a linear 
dependence on a/L. In this way we obtain 

aTpiq) = 0.2374(26) at aspiq) = 0.16535 (4) 

At this point we can see how perturbation the- 
ory can account for this big difference: we know 
from eq. (1) and eq. (2) that the perturbative re- 
lation between asp and axp has a rather large 
1-loop coefficient: 



aTp = asp + 1.5017 x (asp) 



(5) 



It is clear that a naive application of this trun- 
cated expansion could give a not too satisfac- 
tory results for the range of coupling we have 
measured. In fact plugging in eq. (5) the value 
asp = 0.16535 we get for axp the value 0.2064, 
instead of the simulation result given in eq. (4). 

However from previous analyses we know that 
both couplings follow the R.G. flow with the per- 
turbative two-loops /3-function in most of the 
range of scales condidered: this observation can 
easily bring to the conclusion that a perturbative 
relation between aTp{q) and aspifq) (note the 
shifted scale) may give better results provided the 
scale factor / is choosen with care. An intuitive 
choice for / is the ratio of A-parameters, 



/ = Asp/Atp = 0.27620(2) 



(6) 



because with this choice the perturbative relation 
reads 



(XTpiq) = oisp{fq) + 0{{asp{fq)f) 



(7) 



If we now consider the scale qi such that 
«5f(9i) = 0.16535, we get, interpolating our data 
using the effective 3-loops /3-function, 



ocsFifqi) = 0.2289(6) 



(8) 



which is indeed very close to the measured value 
of axp given in eq. (4). By applying this shift 
of the scale to all the measured points, we obtain 
figure 1, in which we can see the two couplings 
nicely superimposing each-other. 
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Figure 1. R.G. flow of the two couplings: scale 
oiasF (diamonds) is shifted. Dotted curve is the 
two-loops universal evolution. 
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